arXiv: 1503.07616vl [math.FA] 26 Mar 2015 


On the determination of a function from its cone transform with 

fixed central axis 


Sunghwan Moon* 


Department of Mathematical Sciences, 

Ulsan National Institute of Science and Technology, 
Ulsan 689-798, Republic of Korea 


Abstract 

A Radon-type transform called a cone transform that assigns to a given function its integral over 
various sets of the cones has arisen in the last decade in the context of the study of Compton cameras 
used in Single Photon Emission Computed Tomography. Here, we study the cone transform for which the 
central axis of the cones of integration is fixed. We present many of its properties, such as two inversion 
formulas, a stability estimate, and uniqueness and reconstruction for a local data problem. 


1 Introduction 


Single Photon Emission Computed Tomography (SPECT), a useful medical diagnostic tool, inspects 
internal organs and produces pictures of their internal processing using the distribution of an isotope. SPECT 
typically provides the information as cross-sectional slice but it is easy to reformat or manipulate into other 
types of images. To obtain the image in SPECT, a gamma-emitting radioisotope is injected into the patient, 
usually, via the bloodstream. This radioisotope passes through the body and is detected by the scan. 

When the conventional gamma camera is used in SPECT, it is common to eliminate large parts of the signal 
since it is mechanically collimated, which means counts only the particles coming from a narrow cone of 
directions and discards the rest. Because of such low efficiency, a Compton Camera was introduced for 
use in SPECT. A Compton camera has very high sensitivity and flexibility of geometrical design, so it has 
attracted a lot of interest in many areas, including monitoring nuclear power plants and in astronomy. 

A typical Compton camera consists of two planar detectors: a scatter detector and an absorption detector, 
positioned one behind the other. A photon emitted in the direction of the camera undergoes Compton 
scattering in the scatter detector positioned ahead, and is absorbed in the absorption detector (see Figure 1). 
In each detector, the position of the hit and energy of the photon are measured. A difference vector between 
two device positions determines the central axis of a cone. The scattering angle ip from the central axis can 
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Figure 1: Schematic representation of a Compton camera 


be computed from the measured energies and electron mass as follows: 


cosip = 1 — 


me 2 AE 
(E — AE)E' 


where m is the mass of the electron, c is the speed of light, E is the initial gamma ray energy, and A E is 
the energy transferred to the electron in the scattering process [1, 14]. Therefore, we get the integral of the 
distribution of the radiation source over cones with central axis /3, vertex u at the position of the scatter 
detector and scattering angle ip. We reserve the name cone transform for the surface integrals of a source 
distribution over a family of cones. 


Many inversion formulas for various types of cone transforms were derived [4, 11, 10, 18, 14]. In particular, 
the cone transform with a fixed central axis was studied in [5, 12, 13, 17, 20]. Cree and Bones derived the 
inversion formula for the cone transform in [5] and Nguyen, Truong, and Grangeat obtained another inversion 
formula in [17]. Haltmeier first defined an ro-dimensional cone transform and found the inversion formula 
in [12]. Jung and Moon discussed a relation between two existing formulas: one derived from Cree and 
Bones and the other from Nguyen, Truong, and Grangeat. In the same paper, they also obtained stability 
estimates for a more general form of the 3-dimensional cone transform in [13]. In this article, we also deal 
this ?r-dimensional cone transform. 


A 2-dimensional cone transform becomes a 17-line Radon transform, which integrates a function along coupled 
rays with a common vertex. This 17-line Radon transform has been studied in the context of Single Scattering 
Optical Tomography [6, 7, 8]. Many works [1, 2, 3, 15, 19] derived inversion formulas for various versions of 
the 17-line Radon transform. 


The definition of the cone transform is formulated precisely in section 2. Section 3 is devoted to elementary 
properties of the cone transform including an analog of the Fourier slice theorem. Two inversion formulas 
are presented in section 4. We describe the range of the cone transform in odd dimensions in section 5. In 
section 6, we show that taking a linear operator on the cone transform is an isometry and discuss a stability 
estimate. In section 7, uniqueness and reconstruction for a partial data problem are studied. 
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2 Definition 


Let f be a function on R 3 with compact support in the upper half space R 2 x [0,oo). We define the cone 
transform by 

27T oo 


Cf(u, 


s) := J J f(u + zs6, z)zd6dz, 


0 0 


for (u, s) = (ui,ii 2 ,s) el 2 x [0,oo) (see Figure 2 (a)). Here 9 = (cos0,sin0) £ S 1 . Let the function / on 
R 4 satisfy /(x,y) = f(x, |y|) for (x,y) = (x\,X 2 ,yi,y 2 ) € R 2 x R 2 and the cone transform of / be defined 
by 

Cf{ u,v) = f /(u+|v|y,y)dy. 

R 2 


Then we have Cf( u, v) = Cf(u, |v|). In fact, making a change of the variables gives 


2n oo 

C/(u,v) = J /(u+ |v|y,yi,y 2 )dy = J J f(u + |v|r0, r)rdrd6 

R 2 0 0 

= Cf(u, |v|). 


Let us consider a natural n-dimensional analog of the cone transform for the function f on R” with compact 
support in the upper half space R ra_1 x [0,oo). As in the 3-dimensional case, the cone transform is defined 

by 


Cf(u, s) := < 


J J f(u + zs6.z)z n 2 dzdS(9), if n > 3, 


S "- 2 o 

OO 


J f(u + zs, z) + f(u — zs, z)dz, if n = 2, 


for (u, s) = (ui, U 2 , • ■ ■ ,u n -i,s) € R n_1 x [0, oo). Here dS(6) is the standard measure on the unit sphere 
S n ~ 2 . When n = 2, Cf is the V-line Radon transform whose integral domain is the set of R-shape lines (see 
Figure 2 (b)). Similar to the definition of Cf for the 3-dimensional case, we define Cf for a function / on 
R 2 ( n " 1 ) with /(x,y) = f(x, |y|) for (x,y) £ R" _1 x R ra_1 by 


C/(u,v)= J /(u+|v|y,y)dy, for(u,v)£ 


R"- 1 


(1) 


Again, we have Cf( u, v) = Cf(u, |v|). Our goals are to reconstruct / (or f) from Cf (or Cf) and to study 
properties of this cone transform. 

Remark 1. The above definition (1) includes the n = 2 case. When n = 2, Cf becomes an integral of / 
along the line perpendicular to _ 

(i, —u)/-\A + H 2 

with signed distance 

u/V^T\ v\ 2 (see Figure 2 (b)). 

In this case the measure for the line becomes 


's/I + \v\ 2 dy. 
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Figure 2: Integration domain (a) a cone (b) a V-line and a line 


Hence we have a relation between Cf and Rf: 


Cf(u, v) = JT+WRft 

\V l + \ v \ 


y/T+W 


where Rf is the regular Radon transform defined by 

Rf(u3,t) = J f{tu: + su: J -)ds, for (w, t) G S 1 x 


Notice that 


Rf 


(!»-«) 


= Rf 


(M) 


yiTKf’ y/T+W Iv^TFF’ y/T+W 


( 2 ) 


3 Elementary properties 

Let <S(R") be the Schwartz class of infinitely differentiable functions f with sup xgR „ |x“c?^f(x)| < oo for any 
multiindices a and /3. We introduce 

5 r (R n_1 x r _1 ) = {/ G ^(IR 2 ^- 1 )) :/(x, y) = /(x, Uy) for any (x, y) G IT -1 x R n_1 

and for all orthonormal transformations U}. 
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Theorem 1. For f £ S r (M. n 1 x R" 1 ) ) we have 

Fi(cm,v) = ^/(€,|v|0 =-F/(£,v|£|), 

where F\(Cf) and Ff are the n— 1-dimensional and 2(n — 1)-dimensional Fourier transforms of Cf and f 
with respect to u and (x, y), respectively. 


Proof. Taking the n — 1-dimensional Fourier transform of Cf( u, v) with respect to u yields 

MCf)(M= J -Fi/(l,y)e i|v|y '«dy = ^/(l,-|v|0, 

R"- 1 

where F\f is the n — 1-dimensional Fourier transform of / with respect to x. Since the Fourier transform of 
a radial function is also radial, we have the assertion. □ 


Note that if / is a radial function on R” and /(x) = f(|x|) , then 

■77(0 = ( 27r)t|^|^H^f(|||), 

where Ff is the n-dimensional Fourier transform of / and 


(3) 


OO 

Hn^f(p)= J f (s)s%jn- 2 {sp)ds, 


where Jk is the Bessel function of the first kind of order k. 

Corollary 1. Let /(x,y) £ iS r (R" _1 x R re_1 ) and f be a function on R" -1 x [0, oo) with f(x, |y|) = /(x, y). 
Then we have 

Fi(cm,s) = (27r)^|^| 2 i 2 S 2 i i H^J- 1 f(C, S |C|). (4) 


Remark 2. When n = 3, (4) was already derived in [5]. 

Proposition 1. The cone transforms C and C are self-adjoint, in the sense that for f,g £ S r 
with /(x, y) = f(x, |y|) and g( u, v) = g(u, |v|), 

J J Cf(u,v)g(u,v)dudv= J J /(x, y)Cg(x, y)dxdy 

R n —1 R n —1 R "- 1 *"- 1 

and 

OO OO 

J J Ci(u,s)g(u,s)s n ~ 2 dsdu = J J f (x, z)Cg(x, z)z n ~ 2 dzdx. 


pn—l w Ttpn-l'v 


(5) 


( 6 ) 


^-l o 


*-i o 


Proof. We start out from 

[ [ Cf{u,v)g(u,v)dudv 


/(u + |v|y, y)dy g(u, v)dudv 


' l ~ 1 R n-1 


v -1 R n-1 R n-1 


J J /(x,y) J 5 (x-y|v|,v)dvdydx 

[“- 1 !”- 1 R ”" 1 

/ / ^ X ’/ 3( x + yl v lw)dvdydx, 


R n — 1 


R "- 1 
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where in the last line, we changed the variables y —> —y. Since / is a radial function in y, we have 


/ / C '^^ U,V ^( U,V 

R n —i R»-l 


)dudv = 


/(x,y) / 

R n_1 


5 (x + y|v|,v)dvdydx, 


which is our assertion. 

Proposition 2. For f,g € S r (M. n 1 x R ra 1 ), we have 

C(f *g)=Cf * Cg. 


□ 


This theorem follows from Theorem 1 and JF(f * g) = (J-f)(J 7 g). 

By the definition of Cf, we notice that Cf is the integral of / over an n — 1-dimensional plane. It is a natural 
idea that Cf can be converted to the 2(n — l)-dimensional regular Radon transform by an n — 2-dimensional 
integration. (Indeed, we showed that the 2-dimensional cone transform is converted to the 2-dimensional 
regular Radon transform in Remark 1.) Let the regular Radon transform be defined by 

Rf(u,t) = J f(tuj + r)dT for (u>, t) G S ”" 1 x M. 

W- 1 - 

Then this can be represented by 

= y/i + KMI 2 J f + t V 1 + dT ( 7 ) 

i “ _1 

for wi 7 ^ 0. Here cj = (cvi,cv 2 , ■ ■ ■ , u> n -i) = (wi, tv') £ S '” -1 and r £ R™ -1 . 


Now we convert the 4-dimensional cone transform Cf to the 4-dimensional regular Radon transform. For 
(a, b) £ R 2 , we integrate Cf(au 2 + b , u 2l v) with respect to u 2 : 

/ Cf(au 2 + b,u 2 ,v)du 2 = / / f(au 2 + b+\v\y 1 ,U2 + \v\y 2 ,y)dydu2 

R R R 2 

= J J f(au 2 + |v|yi - a|v|y 2 + b, u 2 , y)du 2 dy, 


which is equivalent to 


J Cf(au 2 + 6, u 2 , v)du 2 = ((1 + a 2 )(l + |v| 2 )) 1/2 Rf 


(~l,a, |v|, —ct|vl) _ b _\ 

v/(l + a 2 )(l+|v| 2 )V(l + « 2 )(l + M 2 )/' 


(see (7)). Since /(x,y) is radial in y, we have for any a = (cei,a 2 ) £ R 2 , 

f{au 2 + aiyi + a 2 y 2 + b,u 2 ,y) = f(au 2 + agy\ + a 2 y 2 + b,u 2 ,Unyi + Ui 2 y 2 ,U 2 iyi + U 22 y 2 ), 


( 8 ) 
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where U = is any 2x2 orthogonal matrix. Then, we have 

/ / f(au 2 + o.\y\ + a 2 y 2 + b, u 2 , y)du 2 dy 

R 2 R 

= f{au 2 + aiyi + a- 2 y 2 + b,u 2 ,Unyi + Ui 2 y 2l U 2 iyi+U 22 y 2 )du 2 dy 

R 2 R 

= // f(au 2 + (aiUn+a 2 Ui 2 )y 1 +(a 1 U 21 +a 2 U 22 )y 2 +b,u 2 ,y)du 2 dy, 


( 9 ) 


where we changed the variables Uy —> y. Combining two equations ( 8 ) and (9), we have the following 
relation: for \a\ 2 = |v| 2 (l +a 2 ), 


J Cf(au 2 + b,u 2 ,v)du 2 

R 

=«i+« 2 )(i+ M 2 )r i,2 Rf ( 

y 1 /( 1 + o 2 )(i + |v| 2 ) 


v / (i + « 2 )(i + 


By a similar argument, we have the following theorem: 

Proposition 3. Let f £ iSV(R™ 1 x R n 1 ) and (a, a, 6) £ R” 1 x R ra 2 x R. Then we have for |a| = 
|v| 2 (l + |a| 2 ) and u' = (u 2 , U 3 , ■ ■ ■ , u n - 1 ) G R” -2 , 


J Cf(a ■ u' + 6, u', v)du' 


R"- 2 


= ((l + |a| 2 )(l + |v| 2 ))- 1 / 2 J R/ 


(—1, a, a) 


( 10 ) 


V(1 + M 2 )(l + |v| 2 ) y/(l + |a| 2 )(l + |v| 2 ) / 


Remark 3. Taking an inversion formula for the regular Radon transform on (10), / can be recovered from 

Cf. ♦ 


4 Inversion formulas 


Although we already show how to recover / from Cf in Remark 3, we present two explicit inversion formulas 
for the cone transform in this section. 

For k < n we define the linear operator I k by 

ni k m) = i€r fc J7(*). 

The linear operator I k is called the Riesz potential. For / £ S(R n ), T(I k f) £ L 1 (R"), hence I k f makes 
sense and l~ k l k f = f. When I k is applied to functions on M " -1 x R” -1 , it acts on the first n— 1-dinrensional 
variable u or x. Also, for / £ ^(R " -1 x R n_1 ), we have - 7 r i(C'/)(£, v) = .A/(£,£|v|) by Theorem 1, so for 
a fixed v, .Fi(£?/)(£, v) £ <S(R” -1 ) and therefore, I k (Cf) makes sense. 

Theorem 2. Let f £ SV(R ra-1 x R” _1 ) with /(x, y) = f(x, |y|). For k < n — 1, we have 

f = (27 T) l - n If k CI k+1 ~ n F, F = Cf , (11) 
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and 

f = ( 27 r) 1 - n /f fe C/ 1 fc+ 1 _n F, F = C f. 


Proof. We start out from the Fourier inversion formula 


J 1 fc /(x,y) = (2 7 r) 2 ( 1 -") J J | $\~ k Tf(£,r 1 )e i *te i v r 'dtdr, 

JJti-i fn-l 

= (2 7 r) 2 ( 1 - n ) J J \£\ n - 1 - k Ff(£,v\€\)e ix -te iy - v Wd€d-v, 

R "- 1 R 77-1 

where in the last line, we changed the variables r/ —> v|£|. Since Ff(€, rj) and if /(x, y) are radial on 77 and 
y, we obtain 


7 1 fc /(x,y) = (2 7 r) 2 ^ 
which is equivalent to 




J 1 fc /(x,y) = (2^) 2 ( 1 -") J j |Cr- 1 - fc jr 1 (C/)(^,v)e i ( x +lyl v >«^d V . 

R "- 1 R "- 1 


Here we used Theorem 1. Since the inner integral can be represented by the Riesz potential, we have 

/i fc /(x, y) = (2^) 2 ( 1 -") J (I 1 fc+1 -"C'/)(x + |y|v, v)dv 

R 71 - 1 

= (2^) 2 ( 1 -")C(/ 1 fc+1 -"C'/)(x,y). 

and the inversion formula for C follows by applying If k . □ 

Remark 4. Putting k = 0 in (11) yields 

/ = (2ir) 1 - n Cll~ n F, F = Cf. (12) 

When n is odd, equation (12) is actually equivalent to 

/(x,y) = ^-"(-l)^ f A^F(x+|y|v,v)dv, 

R 7 *- 1 

where A u is the Laplacian operator with respect to u. Thus the problem of reconstructing a function from 
its integrals over cones is local in odd dimensions, in the sense that computing the function at a point (x, y) 
needs the integrals over cones passing through neighborhood of that point (x, y). On the other hand, when 
n is even, the inversion (12) is nonlocal because the fractional Laplacian is nonlocal. ❖ 


Combining Proposition 1 and Theorem 2, we have an analog of the Plancherel formula. 

Proposition 4 . Let f,g G S r ( R” 1 x K" x ) satisfy /(x,y) = f(x, |y|) and s(x,y) = g(x, |y|). For any 
k < n — 1, we have 

J J /(x,y)ff(x,y)dxdy = ( 27 r) 1_n J J If k Cf(u 1 v)lf +1 ~ n Cg{u,v)dudv 

Rn-iRn-i R 77 1 R 77 - 1 
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and 


OO OO 

J J i(x, z)g(x, z)z n ~ 2 dxdz = J J I^ k Ci(u, s)/^" +1_rl Cg(u, s)s n ~ 2 duds. 


o f»-i 


0 I »-1 


Proof. By Theorem 2 and the Plancherel formula, we have 

J J /(x,y)s(x,y)dxdy = ( 27 r) 1_n y J /(x, y)Cll~ n Cg(x, y)dxdy 
r -' r - 1 it "- 1 *”- 1 

= (27r) 1_n y y Cf(x 1 y)ll~ n Cg{x,y)dxdy 

i"- 1 1"- 1 

= (2t r) 2 ~ 2n j j |£| fc .Fi ((?/)(£, y)\^\~ k J : i(ll~ n Cg)(^,y)d^dy. 


R n_1 R n " 


Here in the second line, we used Proposition 1. 


□ 


A completely different inversion formula for the cone transform is derived by expanding T\f in spherical 
harmonics, 

OO N(n—l,l) 

Fif(e<p,y) = Yl E ( :p if)ki(Q,y)yik(p), 

1=0 k -0 

and 

OO N(n—l,l) 

Ff(e<p,ii) = E E (FfhifavWkiv), 

1=0 k =0 

where Yik(u>) for uj E S n ~ 2 are spherical harmonics and 

N(n -1,/)= ( ^ + n ~ 3)( oE~ 4) , N(n -1,0) = 1. 


Notice that 


Z!(n-3)! 

y {J r if)ki(Q,y)e~ zy ' ri dy = {Ff) k i{g, rj). 


From Theorem 1, we have the following relation between (J-f)u and ( T\F)ki , where F = Cf: 

W)fci(ftv?) = (V)i-ifev). (13) 


Taking the inverse Fourier transform with respect to v, we have the following theorem: 
Theorem 3. Let f G iS r (R n_1 x ]R n_1 ). If F = Cf, then we have 

n-l 

(Fif)ki(g,y ) = ^ n _i (^)fci(g,yg)- 
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Proof. Taking the inverse Fourier transform of (iFf)ki(Q, vg) gives 


{Fif)ki(Q,y) 


TL — 1 /* Th — 1 /* 

( 2 g 7 r )n-i J (J r f)ki(e,^Q)e ivey dv = (2 g ?r)n _ 1 J (J\ F) kl ( B ,v)e iv e- y dv 

R n -1 r-1 

o n_1 

(27r)"- l (JrF)fc; ( g,yg )' 


where we used (13) in the second line. 


( 14 ) 


□ 


Because of (3), (14) follows 

n+l 3—n 

(TJ) k i(Q,s) = g 2 'Hn-3(TiF));i(g, sg), F=Cf. (15) 

(27r) 2 2 

Remark 5. When n = 3, (15) is already derived in [13, 17]. ❖ 


5 The range 


In this section, we shall determine the range of the cone transform. 

Theorem 4. The cone transform C is a injection o/5,-(R" _1 x R n_1 ) into <S r , c (R n_1 x R ra_1 ), where 
iS riC (R n_1 x R”^ 1 ) = |fe iS r (R n_1 x R n_1 ) : f F(u, v)du is a constant for any v £ R" -1 

l R”- 1 



Proof. Putting £ = 0 in Theorem 1, we have that iff(0,0) = f R „-i Cf(u,v)du is a constant. Hence it is 
enough to show Cf £ 5 r (R" _1 xR” _1 ). To demonstrate this, we show that J r i(C'/)($, v) £ 5 r (R” _1 x R"” 1 ). 
Also, Fi(Cf)(£, v) is infinitely differentiable since F\ (C/)(£, v) = Ff(£, |v|£) = .F/(£, v|£|) and we have 
for any multiindexes a\,a2,(3\, and /? 2 , 

sup{||“ 1 v a2 a| 1 ^ 2 J r i(C'/)(^,v)| : (£, v) £ R”- 1 x R ra_1 } 

= sup{|r i v Q2 ||| l/32| |v|l^l9| 1 ^ 2 J : -/(|,v|||)|($,v) e R" _1 x R"" 1 } < oo. 

Hence, Cf belongs to 5 r>c (R ri - 1 x R" _1 ). □ 


i_ n — 1 n ~ 1 

When n is odd, if " is equal to (— 1)~^~ A u 2 and in this case, we can say more. 

Theorem 5. When n is odd, the cone transform C is a bijection of iS rj o(R ra_1 x R" -1 ), where 


-V, ol 


pn -1 


pn —1 


) = l fes r 


pn —1 


r- 1 ): J /(x,y)dx 


0 for any ygR" 
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Proof. By Theorem 2, we have 


■7 7 1 (C7)(£,v) = .F/(£,v|£|)= f ^if(ty)e~ iy - v ^dy, 

R "- 1 


so J r i(Cf)(0, v) is equal to zero. Therefore, the range of S r r)(&. n 1 x R ra A under the cone transform C is 
a subset of 5 ri0 (K n_1 x R ra_1 ). 


To show C is onto, let F G <S r .o(R n 1 x R™ 1 ). In view of Theorem 2, it appears natural to define 

/ = (27r) 1 - n (-l)^CA^ i F. 

We know that A u 2 F G <S(R" _1 x R" -1 )- By Theorem 4, / belongs to «S r)C (R n-1 x R ra_1 ). In particular, 
•?i/($,y) is equal to 


so for any y G R" - 1 , 


(27r) 1 -"|€r- 1 ^«,y|€l), 


-Ai/(0 ,y)= / /(x,y)dx = 0. 


R"- 1 


□ 


6 An isometry property and Sobolev space estimates 


In this section, we show that I ~ C is an isometry and that the problem of reconstructing from the cone 
transform is well-posed in the following sense : if / satisfying Cf = F is uniquely determined for any F 
belonging to a certain space, the function / depends continuously on F. 

Let L 2 (R 2 ( n- 1 i) be the regular L 2 space. For any 7 G R, let the space L 2 (R " _1 x R" -1 ) and H 7 L 2 { R ™ _1 x 
Rn-i) k e defined by 

L 2 (R n_1 x R”- 1 ) = {/ G L 2 (R 2(n " 1) ) :/(x, y) = /(x, Uy) for any (x, y) G R n_1 x R " _1 

and for all orthonormal transformations U }, 

and 

H^LliW 1 - 1 x R”- 1 ) = {/ G L 2 (R n_1 x R”" 1 ) : ||/|| 7 < 00 }, 

where 

M/Ht = / / l^/(^»7)| 2 (l + l€| 2 Wr7. 

Rn-i Rn-1 

Notice that H 1 L 2 ( R ” _1 x R n_1 ) is the Hilbert space with the norm || • || 7 and H°L 2 ( R "” 1 x R ra_1 ) = 
T 2 (R ” _1 x R"- 1 ). 

Theorem 6. The mapping f —> 2 Cf is an isometry o/LT 7 T 2 (R " _1 x R" -1 ) onto itself. 
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Proof. We start with ||/|| 2 : 


11/11$ = / / l^/(^v)| 2 (l + l^| 2 Wv 

l”- 1 R"- 1 

= / / l^i(C/)(^v/|C|)| 2 (l + |^| 2 )Wv 

p- 1 R"- 1 (16) 

= / / i^i(c , /)«,v)i 2 (i+i€i 2 ri€r- i d€dv 

I”- 1 R n_1 

= (2^) i -"n/r^c/n 2 . 

Here in the second and third lines, we used Theorem 1 and changed variables v/|£| — > v, respectively. It 
remains to prove that the mapping is surjective. It is enough to show that if g £ H 7 L^(R n ~ 1 x R" _1 ) 
satisfies 

J j + \tfyd£dv = 0 

l”- 1 R"- 1 

for all / £ <Sv(R" _1 x R ra_1 ), then 5 = 0. Theorem 1 gives us 

0 = j J ^i5(ev)|^|^7-i(C/)(C,v)(l + |^| 2 )^dv 

l”-i R”- 1 

= IJ ^i5(ev)|^|^^/(C,v|€l)(l + IC| 2 Wv 

l”- 1 R"- 1 

= J J ^i5(€,v/|C|)|CI^^/(^v)(l + |^| 2 )Wv. 

I"- 1 R 1 *- 1 

Since Tf £ i5y(R ra_1 x R” _1 ), Tig is equal to zero almost everywhere, so is g. □ 

Notice that if / £ L 2 (R n_1 x R™" 1 ) with /( u,v) = f(u, |v|), then 

OO 

ll/llo = J J |/(u,v)| 2 dudv = IS" 1 ^ 2 ! J J |f(u, z)\ 2 z n ~ 2 dudz. 

R"- 1 R”- 1 0 p-i 

Let us define iL 7 L 2 _ 2 (R” _1 x [0, oo)) by 

ff 7 L 2 _ 2 (R n - 1 x [0, oo)) = {f : f(u,|v|) = /( u,v), / £ L 2 (R n_1 x R” _1 ), and ||f|| 7 , n _ 2 < oo}, 

where 

OO 

l|f||$,„-2= / [\Fif(t,z)\ 2 m 2 + iyM n ~ 2 dzdz. 

R"- 1 0 

Then if 7 L 2 _ 2 (R" _1 x [0, oo)) is the Hilbert space with the norm || • || 7 i n_2- 

n—1 

Corollary 2. The mapping f I 1 2 Cf is an isometry o/IL 7 L 2 _ 2 (R" _1 x [0,oo)) onto itself. 

The next corollary shows the Sobolev estimates. 
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Figure 3: The upper half unit circle and C-shape lines 


Corollary 3. For each 7 £ R, there is a constant c > 0 such that for f £ S r (M. n 1 x R 71 1 ) with /(x, y) = 

f ( x Jy|), 

||/|| 7 <(2tt)^||C/|| 7+i ^ and ||f|| 7 , n _ a < (27r)^||Cf|| 7+ 


This corollary follows from (16). 

Remark 6. When n = 3, ||f|| 7) i < (27r)~ 1 ||Cf|| 7+ i i i was already discussed in [13]. 

7 The partial data 


From the inversion formula in Theorem 2, / £ <S r (R " -1 x R 71-1 ) is uniquely determined by Cf. However, in 
many practical situations, we know only partial data, i.e., the values of Cf only on a subset of its domain. 
The question arises if this partial data still determines / uniquely. 

Theorem 7. Let f £ 5 r (R n_1 x R" -1 ). The cone transform C/(u, v) is equal to zero for |u| > \/l + |v| 2 
if and only if f(x, y) is equal to zero for |x| 2 + |y| 2 > 1. 

Also, for f (x, |y|) = /(x, y), Cf (u, s) is equal to zero for |u| > \/l + s 2 if and only */f(x, z) is equal to zero 
for |x| 2 + z 2 > 1 (see Figure 3). 


Proof. It is clear that if /(x, y) is equal to zero for any |x| 2 + |y | 2 > 1, then C/(u,v) is equal to zero 
for any v £ R n_1 and |u| > i/l + |v| 2 . From Proposition 3, we know that for |a| = |a| 2 (l + |v| 2 ) and 
(a, a, 6 ) £ R "" 1 x R 71 " 2 x R, 


J C/(a • u' + b, u', v)du' 


= ((1 + N 2 )(1 + | v | 2 ))- 1 / 2 R / 


(—1, a, a) 


\/(l + l a | 2 )(l + | v | 2 ) \/(l + l a | 2 )(l + |v| 2 ) 
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By assumption, we have that 

(—l,a, a.) b 

V(l + |a| 2 )(l + |vP)V(l + |a| 2 )(l + |vP) 

= i/(l + |a| 2 ) (1 + |v | 2 ) f C/(a ■ u' + b, u', v)du' 

R n-2 

is equal to zero for any b/ y /1 + |a| 2 > a / 1 + M 2 > Le -> 

b 

V(i + |a| 2 )(fTM ^) >1 

because the hyperplane {(ui, u') eKx R ” -2 : u\ = a - u' + b} on R " _1 has the normal vector (—1, a) and the 
distance b/y/l + |a| 2 . The support theorem [16, Theorem 3.2 in Chapter II] for the regular Radon transform 
completes our proof. □ 



We study the reconstruction problem for the limited data. The data C/(u,v) is known only for v £ M " -1 
with 0 < a < |v| < b < oo. In order to compute a limited reconstruction, we have to deal with the limited 
cone transform C( aif) )/(u, v) = X a <|v|< 6 ( v )C/( u ^ v )- We define the projection operators by 

P(a,b)f{x, y) = -T r_1 (x a |€l<|r7l<fa|€l ( r 7)- 7r /(€, » 7 ))(x, y) 

and 

P(a,&)f(x,0 = kl^Hn -2 J-lf(|,p))(x, 2 ). 

Theorem 8. Let f £ <S r (R ™ -1 x I"' 1 ). Then we have for k < n — 1, 

P{a,b)Ixf = ( 27 r) 1 -"C'/ 1 fc+ 1 -"C (a , fe) / 

and 

P (Q , b) / 1 fc f = (27r) 1 -"C^' +1 - n C (a , fa) f. 


Proof. By Theorem 1, we have Pi(Cf)(£, v) = Tf(£, v|£), so 

P(P(a, b )I?mtv\Z\) = |l|- fc Xa<|v|<5(v)J-/(|,v|||) = |£|~ k Xa< | v |< 6 (v)Ji {C f) (£, v) 
= \t\- k Fi(c ( a, b) m,v). 

Similar to the proof of Theorem 2, we have 

P (ai6) /f/(x, y) = (2n) 2 ^ J J P(P ia , b) I k m,v)e- i ^ Tl) - (x ’ y) dr 1 d^ 


i-l lttn-1 


= ( 27r ) 


2(1—n) 


?(P(a, b )I k m,v |£|)e- i «- v| « l) - (x -* ) |*| n - 1 dvd£ 


R "- 1 K "" 1 


(27r) 2 ( 1_n ) J J J-i(C' (a , fc) /)(|,v)e-^’l v l«Hx.y)| 4r -i-fc dW|) 


R n_1 R"- 1 


(17) 


where in the second line, we changed the variables r) —> v|£| and in the third line, we used (17) and the fact 
that T\(C( at b)f){€,v) is radial in v. Therefore, we have 

P(a,b)I k f(x, y) = (2^) (1 -" ) J /i fe - n+1 (C (a , 6) /)(x+|v|y,v)dv. 


□ 
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8 Conclusion 


Several types of cone transforms have been studied since the Compton camera was introduced. Here we 
study the n-dimensional cone transform. Two inversion formulas, range conditions, Sobolev space estimates, 
and uniqueness and reconstruction for a limited data problem are presented. 
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